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We assume that the galactic dark matter halo, considered composed of an axionlike par- 
ticles Bose-Einstein condensate tl|, can present topological defects, namely domain walls, 
arising as the dark soliton solution for the Gross-Pitaevskii equation in a self-graviting po- 
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^vq ' tential. We investigate the influence that such substructures would have in the gravitational 

^^1 interactions within a galaxy. We find that, for the simple domain wall model proposed, the 

effects are too small to be identified, either by means of a local measurement of the gradient 
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of the gravitational field or by analysing galaxy rotation curves. In the first case, the gradient 
^N ' of the gravitational field in the vicinity of the domain wall would be 10^'^^ (to/s^)/to. In 

, . ^ the second case, the ratio of the tangential velocity correction of a star due to the presence 

**l^ , of the domain wall to the velocity in the spherical symmetric case would be 10^^. 
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I. INTRODUCTION 

The existence of a mysterious kind of matter, rather different from the usual barionic matter, 
presents itself as a great challenge for modern Physics. This so-called dark matter corresponds to 
almost 23% of the energy density of the Universe ^ and can amount to approximately 90% of the 
total mass in galaxies. 

Recently, it has been proposed that this type of matter can be composed of some kind of 
weakly interacting bosons [3]. When these bosons are spinless they can be identified with axions, 
hypothetical particles proposed in the context of Peccei-Quinn models t4jl. On the other hand, in 
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the case of sub-eV spin-1 particles, they are called hidden bosons or hidden photons [5|, 16(. Axions 
and hidden bosons form a class of particles known as WISPs (Weakly Interacting Slim Particles), 
due to their diminute masses. 

In the last few years, the possibility that the dark matter content of galaxies is in the form of 
a self-graviting Bose-Einstein condensate (BEG) has been considered. Using this approach, the 
authors in [7] were able to relate the mass and the scattering length of an axionlike particle with 
the radius of the galactic dark matter halo. By proposing a new density profile based in the BEG 
features they could construct rotation curves that fit well a sample of galaxies. 

Using the same initial hypothesis, and extending it to spin-1 particles, the authors in [1] showed 
that the mass of the WISP's is constrained by galaxy radii data to the range 10~® — 10~^ eV. 

The next natural step in the identification of the dark matter halo with a BEG is to study the 
possible presence of substructures. BEG's can present a number of different substructures, called 
topological defects, such as vortexes, domain walls, monopoles and textures. 

The existence of these substructures is verified in laboratory experiments, for ultra-cold alkali 
atom gases trapped in a magnetic optical potential, and their features are well studied under these 



circumstances 8h1C1I|. They are observed as a small region (much smaller than the size of the 
condensate) of null mass density in the gas. Mathematically, these defects have origin in zeros of 
the condensate wave-function, stressing the quantum nature of these phenomena in the gas. 



Topological defects have also been studied in the framework of cosmology and gravitation 
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in which they have notable differences from the condensed matter ones. 



or instance, they can 
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have suggested an 



be massive and carry a large amount of energy. Recently, the authors in 

experiment to detect a massive axionic domain wall via magnetic interaction in the context of field 

theory. 

To the authors knowledge, topological defects for a self-graviting BEG have never been proposed 



before. Our goal in the present paper is to explore the possibility that the galactic BEC is endowed 
with a domain wall, a finite region in space where the density vanishes. We are interested in the 
effect that such a structure can have on the galactic dynamics (specially on rotation curves), and if 
it is possible to detect a local domain wall by means of gravitation interaction effects. We restrict 
ourselves to the axionlike particle case. 

This paper is organized as follows. In section [II] we perform the derivation of the density 
functions for the halo endowed with a domain wall, and estimate its width. In section IIIH we 
give an estimate of the gradient of the halo gravitational field in the vicinity of a domain wall. 
In section IIVI we derive a correction term for the rotation curves of stars in spiral galaxies taking 
into account the influence of a domain wall perpendicular to the galactic disk plane (as depicted 
in figured]). Section |V] shows our final remarks. 




FIG. 1. Schematic view of the coordinate system used in this paper. R is the galaxy halo radius. The 
domain wall of width ^ is in a position zq, and is perpendicular to the galactic plane, chosen to be the 
location of the z axis. 



II. DARK SOLITONS 

The zero-temperature mean field energy of a weakly interacting BEC confined in a self-graviting 
potential, V, is given by []| 



E= I d^r 
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where m is the mass of the particle composing the condensate, a is the s-wave scattering length 
and ilj{r) is the condensate wave function, satisfying f d^r\ip\'^ = N with N being the total number 
of particles. 

The mean field dynamics of the system is described by the Gross-Pitaevskii (GP) equation 

in^!t^ = (-|^V2 + y(r) + g|^(r, t^^ ^(r, t) , (2) 

where g = ATrh'^a/m. 

Some of the solutions of equation ([2|) may be quantized vortexes or dark solitons. These func- 
tions are topological defects in scalar BECs, in which the density vanishes due to the topological 
constraint on the phase of the wave function. 

In order to investigate the topological defect that corresponds to the domain wall that can 
appear perpendicularly to the z direction, we coupled the topological defect with the ground state 
of the galactic condensate in the Thomas-Fermi approximation [7| in the form 

tpir, t) = il^TFix, y, z)(l){z, t) , (3) 

where tpTFix^y^z) = tl^TF{r) is the Thomas-Fermi solution for the GP equation (with r^ = x^ -H 

y^ + z'^), 
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for r > R 



with k = ^^/Gm^/Wa, R = vr/fc is the condensate radius and po is the central number density of 
the condensate. (t){z, t) corresponds to the topological defect solution. 

We are interested in characterizing the defect by its position z, hence we eliminate the Thomas- 
Fermi solution in GP equation, as well as its dependence on x and y coordinates, by multiplying 
([2]) by ijj^p and integrating it in these coordinates, obtaining 

^h^-^^= (-^^,-v{z)^ + V + 9{zmz,t)A<P{z,t) , (5) 



dt \ 2m dz"^ dz 



where r/(z) = —^^ ( i+ cos (kz)) ^® ^^ extremely small factor and the term it couples to can be 
neglected. 

The effective interaction parameter, g{z) = gpof{kz)/2, is proportional to the central density 
and the form factor 



where x = z/R. 



"^ 



fi^) 



1 + cosfvrx) 



(6) 




FIG. 2. Form factor of the condensate. 



We assume that the width of the topological defect is much smaller than the size of the con- 
densate. In this situation we can consider the particle number density and the effective interaction 
parameter as almost constants in the defect vicinity. As the self-graviting potential obeys the 
Poisson equation, we can approximate the potential as V{r) ~ Vq and, therefore, we substitute 



(l){z,t)=u{z,t)e-'^'''^'' 
in equation ([5]) to obtain the one-dimensional Gross-Pitaevskii equation for u{z, t) 
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du{z,t) 
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+ g'\u{z,t)\'^]u{z,t) 



dt \ 2m dz"^ 

where g' is supposed be locally constant. Equation ([8|) has the solution 
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(z, t) = e-^/^^/Hanh 
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where fi = g' = g{zQ). This solution is called a planar dark soliton, describing a domain wall at 
z = zq, since the density vanishes at that point. 



TABLE I. Values for the healing length ^ using masses and scattering lengths obtained in 
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a (fm) 
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10-5 
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The quantity S,, called healing length, is related to the width of the domain wall. It is possible 
to show that it is a function of the parameters of the condensate in the form 

^ (10) 



^/47rpoaf{kz) 

The density function for the domain wall pow = l^(-2^)^)P is shown in figure [3l 
1.0 f 
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FIG. 3. Density function for the domain wall located at the origin of the coordinate system and with a 
healing length of 0.1 (in arbitrary units of length). 



The density function 



Pz 



(r,t)| dxdy 



(11) 



for the condensate with a domain wall is depicted in figure HI 

Local dark matter density measurements indicate a value of OAGeV/cm^. Using that informa- 
tion and the masses and scattering lengths for an axionic dark matter halo as suggested in [l|, we 
can infer the order of magnitude for the heahng length of the domain wall. The results are shown 
in table [H 
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FIG. 4. Density function (in units of k"^ /{2pq)) for the condensate endowed with a domain waU. The width 
of the domain waU has been made large in order to facilitate visualization. 

We can see that the healing length decreases very quickly with the particle's mass, becoming 
negligible for larger masses. The value for a mass of 10~^ eV is already beyond any physical 
significance at galactic scales. 

III. GRAVITATIONAL EFFECT IN THE VICINITY OF THE DOMAIN WALL 



We proceed now to the calculation of the effect of a domain wall located on the galaxy disk, 
more specifically crossing Earth's position, on a test body (e. g., a satellite). 

In the presence of a domain wall, the total density distribution is not symmetrical, then the 
gravitational effects are distinct from the case of the halo density without the domain wall. We 
intent to estimate this difference by analysing the movement of a massive test body crossing the 
domain wall. 

The gravitational field on the body is given by the solution of the equation 



V-5 = -47rG^, 
where g is the mass density which can be related to the wave function by 

Q{x,y,z) = m\i;{r,t)f . 



(12) 



(13) 



The gravitational effect in the test body will be maximized if this body is moving along the 
z-axis and between the border and the center of the domain wall. By the symmetry of the spatial 



configuration, the gravitational field has only a z-direction component. In this case, the equation 
to be solved is 

d , . . ^ sm( kz) , ■) ( z — zn\ ,^ ,, 

-,A.) = -A.G,np„^V..Vf [-^) . (14) 

The difference in the gravitational field is given by 

gz{z,) - g,{zo - i/2) ^ -47rGmpoe^^^^^ (^/2 tanh (^^) - l) , (15) 

where xq = zq/R is the domain wall position relative to the galaxy radius. 

For the Sun's relative position, xq ~ 0.5 and the gravitational effect exerted in the massive body 
crossing the topological defect is of the order of 10"^*^ m/ s^ , for a healing length of the order of 
10^ m and the gradient of the gravitational field is 10^'^^ {m/ s^)/m. Because the Earth's movement 
(along with the Sun) in the galaxy has a velocity of ~ 10^ m/s, this effect in the vicinity of our 
planet could only be detected by an experiment with a precision greater than 10"'^^ m, which is 
far beyond present day technological capability. 

IV. TANGENTIAL VELOCITY CORRECTION 

When the domain wall is present, the gravitational field presents tangential components. How- 
ever, the projection of the gravitational field in the tangential direction is much smaller than the 
projection in the radial direction even near the domain wall. Then, we can neglect the tangential 
components and assume that the gravitational field is radial and given by 

^|-r-V(O = -4vrG^(^,0,</>), (16) 

r^ or 



where 



/ ^ ,x sin(A;r) , n /r cos(6') — zn'i /^«n 

Q{r, 0, 4,) = "^Po— ^ tanh2 ( \1 ) . (17) 



Using Gauss theorem in the equation (fTUj) . we obtain 



9ri^) = -2 ' (18) 



where the mass profile of the dark condensate galactic halo is, 



MDM{r)= / ^(r,^,0)d^r, (19) 

Jv 



with V the volume of a sphere with radius r. 



Equation ([18]) allows to represent the tangential velocity v^g{r) = rgr{r) of a test particle moving 
in the halo as 

Vtgir) = v'^^^ir) - v'^„„{r) , (20) 

where 

"Ll--) = ^^ (^ —(*-•)) (21) 

is the squared tangential velocity for the spherically symmetric case (already obtained in [7]) and 

2 47rGmpo f r- ^ cos{kzQ) - cos{kr)\ 
Vcorrir) = p ^V27T-e{r - Zq) j (22) 

is the correction in the squared velocity due to the presence of the domain wall. Q{x) is the 
Heaviside step function. 

f^o^j. is proportional to ^/R <^ 1, then the correction is maximal when the wall is located near 
the center of the galaxy. As the domain wall width is always many orders of magnitude smaller 
than the radius of the halo, this correction is also small. 

In figure [5] both terms of (I20p are shown, to stress the difference in magnitude they present. 
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FIG. 5. Tangential velocities (in km^ / s^) for the BEG dark matter halo (dashed line) and the domain wall 
(solid line). The wall's relative width ^/R has been chosen as 0.05 to allow easy visualization of both curves. 
It is possible to see that the diference between the curves is very large. Here x = r/R. 



With the addition of a barionic matter term (after choosing an appropriated barionic matter 
density profile), equation ()20p can represent a rotation curve for stars in a spiral galaxy. The term 
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FIG. 6. Tangential velocity correction due to the influence of the domain wall with a relative width of 0.01. 
The curves are related to walls located in increasing relative distances xq from the center of the galaxy. 
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v^.g had already been obtained in [7], and it was found to fit observed rotation curves for a number 
of galaxies. Our correction, v^^^^, because of its small magnitude, cannot be detected in these 
types of rotation curves. As the factor ^/i?, for typical galaxies, would amount to about 10~^^, the 
ratio between the correction and the tangential velocity for the spherical symmetric case would be 

VcoTrhss ~ 10"^. 

The small efect that a domain wall would have in the galactic dynamics can be explained when 
we take in consideration the mass that could fill a disk the same size as the domain wall in a galaxy 
similar to the Milky Way {R. w 10 /cpc). This mass would amount roughly to lO^'^ kg^ or about the 
mass of the Moon. 



V. CONCLUSIONS 



By assuming that the dark matter halo in galaxies is composed of a condensate of bosonic 
particles (with axionlike properties), as previous works hypothesized, we were able to model one 
type of substructure in the halo, in the form of a topological defect known as domain wall, derived 
from a dark soliton solution for the Gross-Pitaevskii equation in a self-graviting potential. 

Because other types of topological defects (such as vortexes, monopoles, textures and ring 
solitons) would occupy a smaller volume in the halo, and therefore would have a smaller influence 
on the total dark matter density, we decided to restrict to the study of domain walls. Even in this 
case, the magnitude of the effects are too small to be subject to detection by present methods, at 
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leas, for the choice of para.„e.e. ,.„ain,y the heahn, length t which depend, on the .a. and the 

scattering length of the axionhke dark matter particle estimated in [1] ) and simplifications we have 
made here. For the local gravitational interaction, we have a gradient in the field of the order of 
10^'^^ {m/ s^)/m. The correction factor on the velocity rotation curves for stars in spiral galaxies 
is typically of the order of 10^^. 

However, there may exist some kind of cumulative effect that renders the influence of domain 
walls considerable in a system with a larger number of topological defects or a greater dark matter 
density. They also may be important in other phases of galaxy evolution. 

The main result of this work is the implementation of a methodology for the inclusion of 
topological defects in a quantum gas dark matter halo. The sequence of this study implies, for 



ity of the domain wall, the 

1 

1| ) , and the manifestation 



example, the determination of the dynamical and thermodynamical stabi^ 

extension of the method for a spin-1 particle condensate (as suggested in 

of such defects in a fermionic quantum fluid, among other possibilities. These issues will be the 

subject of future work. 
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